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Abstract 

We propose a new framework for transverse-momentum dependent parton distribution functions, 
based on a generalized conception of gauge invariance which includes into the Wilson lines the 
Pauli term ~ i 7 ^^, 7^]. We discuss the relevance of this nonminimal term for unintegrated 
parton distribution functions, pertaining to spinning particles, and analyze its influence on their 
renormalization-group properties. It is shown that while the Pauli term preserves the probabilistic 
interpretation of twist-two distributions — unpolarized and polarized — it gives rise to additional 
pole contributions to those of twist-three. The anomalous dimension induced this way is a matrix, 
calling for a careful analysis of evolution effects. Moreover, it turns out that the crosstalk between 
the Pauli term and the longitudinal and the transverse parts of the gauge fields, accompanying 
the fermions, induces a constant, but process-dependent, phase which is the same for leading and 
subleading distribution functions. We include Feynman rules for the calculation with gauge links 
containing the Pauli term and comment on the phenomenological implications of our approach. 
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I. INTRODUCTION 



Parton distribution functions (PDF)s are the key nonperturbative ingredients of com- 
pletely inclusive QCD processes, like deeply inelastic scattering (DIS). The process- 
dependent hard-scattering part of such processes can be calculated order by order in QCD 
perturbation theory on account of the hard scale of the process Q 2 ^> Aq CD . Though the 
determination of the initial PDF requires the application of nonperturbative methods, its Q 2 
evolution is controlled by renormalization-group (RG) evolution equations with anomalous 
dimensions calculable within perturbative QCD. 

This simple picture changes significantly when one considers semi-inclusive processes, 
like semi-inclusive DIS (SIDIS), or the Drell-Yan (DY) process in hadronic collisions, in 
which hadrons are detected in the final (initial) state with a sizeable transverse momentum. 
In that case, one needs information about the generation of the transverse momentum P^ 
of the final (initial) hadrons, e.g., by means of the transverse- momentum distribution of 
the partons. This mechanism is believed to be dominant at small P± <C Q, while at large 

~ Q, the transverse momentum P\ is produced by the perturbative gluon exchanges. The 
second mechanism, as well as the relationship between the two in the intermediate region, 
are outside the scope of the present work. In any case, integrated PDFs of leading twist 
are not sufficient to describe semi-inclusive processes. One therefore introduces transverse- 
momentum dependent (TMD) PDFs which keep track of the intrinsic transverse motion 
of the partons inside the hadrons and reveal this way fine details about their substructure 
(pioneering works are QH— see also QJloj] and Refs. cited therein, and for a review). 



The introduction of TMD PDFs, though intuitively clear and physically appealing, still poses 
serious challenges. The first problem is related to the TMD factorization: Its status beyond 



leading twist (and to all orders) is far from being satisfactory at the moment 12H15I]. Next, 



there is a possible non-universality of TMD PDFs entailed by the extremely complicated 
and often process-dependent structure of the gauge links, 1 see Refs. 16l4l8ll . Finally, in the 

nnnrn 

light-cone gauge, extra divergences appear that have to be properly treated |4Ll5Lll9Ll20] — in 
contrast to the integrated case. In the present paper, we focus on the last two issues. 

In the integrated case, the parton density fi/h(x,Q 2 ) describes the probability to find a 



1 These are path-ordered exponentials of the gauge field, needed to render the definition of PDFs gauge 
invariant. 
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parton i with longitudinal momentum fraction xP + inside hadron h with momentum P, 
and can be given a gauge-invariant definition in terms of the gauge link (Wilson line) [see, 



for instance, 



2l|] 



[£-;0-|r]=Pexp 



ig / dz»A*Jz)P 



o-[c] 



(1) 



for a contour C along the light-cone, where the path-ordered exponential A a ^ refers to the 
(gluon) gauge field. The renormalization of the integrated PDF obeys the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) 22|, |23[ evolution equation with its integral kernel 
being related to the anomalous dimension IY^ = S q r g + 2r end , where r en d is the endpoint 
anomalous dimension of the integration contour C in the gauge link (see 24J for a more 
detailed discussion of this issue and |25l-l27l] for the original derivations and earlier references). 



It was pointed out in ua 12a i2si that a completely gauge-invariant definition of the 
TMD PDF in those gauges in which Aj_ does not vanish at infinity has to include also 
transverse gauge links. Hence, in the light-cone gauge A + = 0, applied in conjunction with 
^'-independent pole prescriptions (like the advanced, retarded or principal-value prescrip- 
tion) in order to avoid singularities at q + in the gluon propagator, the transverse gauge links 
receive radiative corrections that can be associated with a cusp-like junction point at light- 



cone infinity 



24 



30 



3l|. The emerging cusp anomalous dimension 32] has to be removed, 



if one aims to recover the results valid in covariant gauges. To this end, a redefinition of 
the TMD PDF was proposed by two of us 



24 



30] which involves a soft factor, termed R, 



consisting of two eikonal lines evaluated along a particular gauge (integration) contour with 
a jackknifed path segment in the transverse direction (see next section). Note that the in- 



troduction of the soft 
instance, Refs. 



17 



33 



'actor can be justified from a different point of view as well — see, for 



34J , where the soft factor was used to take care of rapidity divergences 



in covariant gauges. The anomalous dimension related to the ultraviolet (UV) diver gen ces 



(e.g., pole terms in 1/e in dimensional regularization) of the soft factor was found |24j, |30j to 
exactly compensate (at the one-loop order) the cusp anomalous dimension of the transverse 
gauge link, hence, ensuring the independence of the (redefined) TMD PDF from artificial 
contour-generated anomalous- dimension artifacts. 

More recently 35], we have shown that this factorization scheme remains valid also for 



the Mandelstam-Leibbrandt pole prescription 361. l37j|. which is q -dependent. In that case, 



the UV divergent part of the soft factor reduces to unity, while the transverse gauge link 
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does not give rise to a defect of the anomalous-dimension that has to be compensated. As a 
result, the TMD PDF has the same anomalous dimension as in covariant gauges, rendering 
the proposed definition of the TMD PDF gauge and pole-prescription independent. 

The basic tenet in the gauge-invariant formulation of hadronic quantities, like TMD 
PDFs, is to use a gauge link with an exponent which contains only the gauge field A. 
However, this is only the simplest (or minimal- twist) possibility which pays attention to 
the fact that color vectors cannot be compared at a distance. Because the gauge poten- 
tial is spin-blind, one should actually include into the gauge link an additional term 
proportional to the gluon tensor — called the Pauli term — which can accommodate the 
direct spin-dependent interaction in accordance to the Lorentz group. This term repre- 
sents the minimal coupling of a spinning particle to an external field and may become 
important for non-trivial contours, while additional terms of still higher twist are not pro- 
hibited but are relatively power-suppressed. Thus, the gauge links will be generalized to 



take into account the Pauli contribution 



381 



Fjt»Siu>, where s tw = (V 4 )^^], nor- 



mally ignored. This means that in order to accommodate spin-dependent interactions in a 
manifestly gauge-invariant formalism, one has to include the following path-ordered expo- 
nential Pexp j — ig J °° da S ' llv F£ h '(ua)t a ~\ . The graphic illustration of this concept is depicted 
in simple contextual terms in Fig. [1] which shows a generic process with gauge links that 
contain the Pauli term — codified by small rings around the double lines which stand for the 
conventional gauge links. 

X _ 



(oo-,0±) 



A 



V 



(oo- £l) 



FIG. 1: Graphical representation of a generic TMD PDF in coordinate space with the spin- 
dependent Pauli terms included in the gauge links. The (hidden) contour obstruction at transverse 
and lightlike infinity ("earth" symbol) is represented by x. 
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It is expected that any effects of such spin-dependent terms should be non-vanishing only 
in the case of (at least) transverse gauge links off the light-cone. In the integrated collinear 
PDFs, the Wilson lines are one- dimensional in the sense that the paths of the integration 
reduce to lightlike rays. Going beyond the fully collinear picture, which is unavoidable in 
semi-inclusive processes, one must make use of gauge links which involve more complicated 
integration contours, which have at least one additional — transverse — dimension. Note that 
in the case of integration paths off the light-cone — used to regularize rapidity divergences 
4( — one has, in fact, even more non-trivial contours because they contain, beyond the minus 
light-cone segments, also plus components. Therefore, in the TMD case, effects related to 
the spin transfer from the starting point, say, (0~,OjJ, to the terminating point (£ _ ,£j_) 
may become (at least, in principle) apparent due to the non-trivial structure of the contour. 
Hence the cross-talk between a pair of quantum fields at distant points off the light-cone 
will contain spin-dependent Pauli terms which are of higher twist order with respect to the 
spin-blind ones containing the gauge potential. Our analysis reveals that the inclusion of the 
Pauli term, although non-visible in the completely unpolarized TMD PDFs, can produce 
non- vanishing effects in a number of polarized distributions, in particular, those responsible 
for time-reversal-odd phenomena. 2 

Adopting this encompassing concept of gauge invariance, questions arise whether the 



definition of TMD PDFs, we proposed before in Refs. |24j . |30|, has to be modified and 
whether the inclusion of the Pauli term has phenomenological consequences — as already 
indicated. The first issue is related to the question whether spin-dependent terms affect the 
factorization schemes discussed in our previous works, while the second addresses possible 
changes of the RG, i.e., evolution behavior of TMD PDFs. The present work is devoted to 
the clarification of these issues. 

The rest of the paper is organized as follows. The pivotal Sec. [Ill argues that the correct 
treatment of spin degrees of freedom in the TMD PDFs necessitates the inclusion into 
the gauge links of the Pauli term. This contribution describes the interaction between 
spinning particles and the gauge- field strength and leads to a generalization of Eq. (CQ). Its 
implications are worked out in Sec. Ill II The calculation of virtual gauge-field correlators for 
the leading-twist distributions, as well as for those of subleading twist, is carried out in Sec. 



We thank A. V. Efremov for important comments on this point. 
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HVl whereas those related to fermions are discussed in Sec. [V] Section [VI] is concerned with 
the consideration of contributions to the TMD PDFs stemming from real-gluon emission. 
Finally, in Sec. IVI!} we summarize the results and present our conclusions. To go further 
with QCD calculations with gauge links, which include the Pauli term, we develop a set of 
Feynman rules and display them in Fig. 



II. INCLUSION OF SPIN EFFECTS 

The TMD PDF for an unpolarized/polarized quark of flavor i in an unpolarized/polarized 
target h following our generalized concept of gauge invariance reads 

fU*,k ± ) = ^Tr J dk~ J-0yL^ ik< (h MO 

x [[£- , £l\ oo-, £±\H°o-, £ ± ; oo-, oo ± ]f ■ V ^ 
x [[oo-, ocu; oo", 0j.]][[oo-, 0±, 0", 0^(0)^) -R 

where T denotes one or more 7 matrices in correspondence with the particular distribution 
in question, and the state \h) stands for the appropriate target. In the unpolarized case, 
we have \h) = \h(P)), with P being the momentum of the initial hadron, whereas for a 
(transversely) polarized target the state is \h) = \h(P), S±). The "enhanced" gauge links 
[[£ 2 ;£i]] and the soft factor R will be defined shortly. 

An important comment about definition (|2j) is here in order before we proceed. We started 
from the "fully unintegrated" correlation function, which depends on all four components 
of the parton's momentum 13|, |39[. Thus, the TMD PDF is obtained after performing the 



k integration, which formally renders the coordinate £ + equal to zero: 

dk~e- ik ~t + = 2tt(5(£ + ) . (3) 



However, one must be careful: This operation may produce additional divergences because, 
carrying it out, all quantum fields involved (quarks and gluons) will be defined on the light 
ray £ + = 0. This means that the plus light-cone coordinates of the product of two quantum 
fields will always coincide. To avoid this, we will regularize this singularity in what follows 
by taking into account that a particle, once created at the point £ + = 0, will be reabsorbed 
(destroyed) with the same probability at (potentially very distant) points ± A/2, where 
A ~ l/p~ ~ p + /2M 2 is the uncertainty of determining a point along the plus direction. 



In other words, we have to sum (average) over all indistinguishable possibilities in order 
to get the correct answer in the quantum mechanical sense. For instance, the regularized 
two-gluon correlator is written as 



A/2 



A/2 



dt (^(o+ r,e±M"(*,?- &)>,, 



A— i>oo 



T J (2tt) 4 
x 2nd (q~) D^(q) 



(4) 



whereas without regularization, the corresponding term ~ J dq~ D^ v {q + , q~ , q ± ) would face 
unphysical UV divergences. 

The constant T ~ l/p + (so to say the "length" of the plus ray) will drop out from all 
final results, provided a suitable parametrization of the vectors along the contour integral is 
adopted. This is crucial for the enhanced gauge link which includes the Pauli term, since the 
latter is not reparameterization invariant — in contrast to the usual gauge link. Therefore, 
we make use of the following reparameterization of the (initially dimensionless) constant 
vectors that define the motion along the line integral: 



n n^ u »= P n P 



l± -> P + L 



(5) 



which implies boosts in the collinear directions. Note that the plus-component of the mo- 
mentum, p + , is large in our kinematics and is the only mass scale entering the above repa- 
rameterization. Thus, the uncertainty of determining a position along the plus ray in Eq. 
(jj]) is very large, namely, 

A 1 1 

A ~ — = — . 

u* p 

while it is very small along the minus or the transverse directions: 

1 1 1 



A 



r\j — r\J 

u \l 



We can now define the enhanced lightlike gauge link along the x direction: 



[[<xT, Ox; 0- 0J] = Pexp 



-w 



poo poo 

/ da w M A^{ua)t a - ig da S^F^(ua)t a 
Jo Jo 



■ (6) 



An analogous definition holds for the x + direction by making the replacement u — > u*. On 
the other hand, the enhanced transverse gauge link is given by 



[[oo ,cxdj_;oo , Oj 



Pexp 



poo poo 

-ig / dTl ± -A a ± (lr)t a -ig / dTS^F^(lr)t a 
Jo Jo 



(7) 



S 



(0+ oo- £_) 




(-00+, 0", ± ) 



FIG. 2: Integration contour associated with the soft factor R in Eq. Q. The cusp angle \ is 
explicitly indicated. 

where the two-dimensional vector I = Zj_ drops out from all final results, and the Lorentz 
generators for the spin are defined by S^ v = (1/4) [7^, 7^]. Note that the path ordering, 
denoted by V in the compound expressions above, means 



Tl 



V[. ..] = ! + J dr{P exp dru- Aj gS ■ J*(wri)P exp dru ■ A 

POO PT2 / POO \ / l>T2 

+ / dr 2 / driPexp I / dru ■ A I ■ gS ■ J r (uT2)Vexp I / dru- A) (8) 



n 



x gS ■ J-~(uti) ■ V exp I J dru ■ A J + . . . , 

where we have used the following convenient abbreviations: A = J2 a A a t a , u-A = J^^u^A^ 1 , 
S ■ F = ^ v S^ U F^ U , and the path ordering inside Eq. ([8]) is the usual one. It becomes 
obvious that the enhanced gauge links, defined above, and the standard ones fulfil the same 
gauge transformations. 



The soft factor R in Eq. §8§ — introduced in [24j, |30| with the aim to remove the defect of 
the anomalous dimension of the TMD PDF — may, in principle, be upgraded to include the 
tensor term as well. This amounts to the following expression 



R(p + ,n-\C,Z±)=Tr (0 



Pexp 



ig ds(-A(()+ig dsS-T(() 



C-C.1 



x Pexp 



ig ds(-A{S + 0-ig dsS + 



Cc.1 



(9) 
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where ((s) = d(/ds, V denotes anti-path ordering, and the integration contour C cusp is 



the same as that employed in 24], |30| (see Fig. [2] for an illustration). Note in this context 



that the soft factor was introduced before (without the Pauli term) in Refs. [33|, |34j with the 
purpose to control rapidity divergences of non-lightlike Wilson lines in covariant gauges. The 
soft factors in both approaches are multiplicative renormalization eikonal factors, though in 



33 



34J the contribution from the gauge link at infinity is not considered owing to the use 



of a covariant gauge. 



III. INFLUENCE OF THE PAULI TERM 



Before we focus our attention to the specific implications of the Pauli term, let us first 
summarize the key features of the proposed scheme. The usefulness of Eq. ([2]) derives 
from the fact that by virtue of the soft factor R all gauge-dependent anomalous-dimension 



artifacts, potentially contributing to the TMD PDF, are absent ab initio [24J, [30J so that, 
integrating over the transverse momenta, one obtains a PDF which is controlled by the 
DGLAP evolution equation 20|, |24j with the usual anomalous dimension. Moreover, to this 
definition all pole prescriptions adopted to evaluate the gluon propagator in the light-cone 
gauge are fungible [35]. 

To study the effects of the spin-dependent terms, induced by the inclusion of the Pauli 
contribution, it suffices to take them into account only in the fermionic part of Eq. (|2]), 
leaving the soft factor unmodified. The justification of this treatment is based on the fact 
that the structure of the soft factor is practically prescribed by the RG properties of the 



30 



351 ] . To be more specific, it was 



unsubtracted TMD PDF, as shown in detail in Refs. 
found there that the particular contour C cusp in the soft factor, depicted in Fig. [2] pertains 
to the cusp-like UV singularities of the fermionic part of Eq. ([2]). Another argument of 
retaining the original form of the soft factor unchanged is provided by the requirement that 
it should be boost invariant (see, e.g., Ref. jl3|). Given that the Pauli term is not invariant 
under scale transformations, we refrain from including it into the soft factor in the present 
investigation. From the calculational point of view, the above argument is related to the 
fact that, in the absence of any Lorentz structure, the spin-field interaction cannot produce 
nontrivial results for integration paths without self-intersections — this will be considered 
elsewhere. 
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Carrying out the k integration in Eq. fl2]) and leaving out the soft factor R, one obtains 
the following unsubtracted TMD PDF 



f[ /q (x,k ± ) = ^Tr J ex P {~ ik+ t, +ik±-£±) (p,s|V>i(f >^x5 00 >£J] f 

x [[00", ooJ^nfocT, 00 ± ; oo~, 0_l]][[oo~, 0_l; 0~, Oj] 

xVi(o _ ,o±)b,s> • 

(10) 

For our concrete calculations to follow, we consider matrix elements between quark states 
having momentum p and spin s: \p, s). Using the fermionic density matrix (p = p ■ 7) 

u(k) <g> u(k) = -{% + m) (1 + 75S) , s 2 = -1 , (11) 



with the spin vector = (s + , s , s±) being given by 39] 



\ m ' 2mk + ' m J \ ' /c + ' 



one obtains in the tree-approximation (indicated by the subscript for a°) 

/( r o)0r,A^) = ^Tr[(p + m)(l + 7 5 5) T] 5(p+ - xp+)5^(k ± ) . (13) 

In particular, for the unpolarized TMD PDF with T = 7 + , one has at leading twist two the 
expression 

/f ee /Jt(x, fcj.) = ^Tr [(p + m) (1 + 75s) 7 + ] Kp + ~ xp + )5^(k x ) 

1 (14) 
= 5{l-x)5 (2 \k ± ) . 

On the other hand, the helicity and the transversity distributions are given, respectively, by 

f$ v (x,k x ) = 6(l-x)6®(k ± ) A, 

(15) 

where A denotes the helicity and the transverse spin of the parton quark i. 

To continue this kind of calculation beyond the tree level, we have to expand the product 
of the enhanced gauge links and retain all terms contributing up to 0(g 2 ). That is, employing 
the light-cone gauge A + = (A ■ n~) = 0, we have to evaluate 

[[ocT, 00 j_; ocT, Ox]] • [[oo~, 0_l; 0~, 0_l]] a +=o = 1 - ig (Ui + U 2 + U 3 ) 



-g 2 (U A +U, + ...U W ) 

11 



(16) 



TABLE I: Individual virtual-gluon contributions appearing in the evaluation of the product of the 



gauge links in Eq. (|16p up to the order 0(g 2 ). 



Symbols 


Expressions 


Figure [3] 


Value 


U\ 


Cdr I ■ A(It) 


(a) 


± 0, [241 


Ui 


Jo°dT S • J~(ut) 


(b) 


7^ 0, see text 


U 3 


r n °°dr S ■ T(It) 




0, Eq. ([21]) 


u 4 


r n °°dr (7 (icr (Z • A(It)) (I ■ A(la)) 

JO JO V v J J \ v // 




0. [241 


U 5 


CdrJS'cfcrCl-^Ir)) (5-^(Ia)) 




0, Eq. @U 


U G 


SS°dT£dff(S-niT)) (l-A{la)) 




0, Eq. ([H]) 


U 7 


Io° dT Jo da ( S ' jr (. UT )) ( S ■ Fiua)) 


(c) 


0, see text 


U 8 


$™dT$Zdo(S-F{lT)) (S-F(la)) 




0, Eq. <[21|) 


U 9 


J °°dTf °°da (l-A(lr)) (S-T(ua)) 


(d) 


7^ 0, see text 


Uio 


JS°*rJZ°d* (S-Hlr)) (S-T(ua)) 




0, Eq. (EU 



where the individual contributions entering this equation are compiled in Table |U 3 In these 
expressions we used for the sake of convenience the longitudinal vector = p + n~ which 
has units of mass instead of the dimensionless light-cone vector n~, cf. Eq. (jHJ). 

The entries in Tabled call for some comments and explanations. First, the fermion fields 
in the definition of the TMD PDF given by Eq. fTTUT) are Heisenberg field operators, meaning 
that we have to use 



drj ipAip 



= / d\ ^(vhMv^iv) . 



(17) 



Therefore, the O(g) contributions ^1,^2,^/3 should be contracted with the quark-gluon 
interaction terms J drj ipAijj , originating from the Heisenberg fields ( fTTI) . in order to give 
rise to the one-gluon exchange graphs (a) and (b), which are of 0(g 2 ). Second, all virtual- 
gluon terms iVj ee (Ui) with i = 1, ... 10 produce contributions of the following generic 



The nonlinear part of the gluon tensor does not contribute in the considered order of the coupling. 
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form 4 

~ / ^f4% ex P i-'(P + - k + )C+iks. ■ £l) \tt [(p + m)(l + 75 5) T (iV)] 
= 5(p + - xp + )<5 (2) (fc±) ^Tr [(p + m)(l + 75 s) T (N)] . 
Their Hermitean conjugated (mirror) counterparts contribute terms of the form 

~ 5(p + - xp + )5 {2) (k 1 _) ^Tr [(p + m)(l + 7 5 s)(A^) t T ] . (19) 

The Dirac structure of the quantities (Ui) is nontrivial owing to the spin- dependent terms 
from the gauge links — which we will show explicitly below. On the other hand, the con- 
tributions of the real-gluon exchanges, stemming from the contractions of the gauge fields 
belonging to different planes in the £-space, will be considered in Sec. [VTl 

It is obvious that an analogous expansion has to be carried out in Eq. (fTUl) also for the 
product of the gauge links [[£", oo~, £x]]* • [[oo - , oo~ , ooj_]]t. Let us emphasize that 
the various contributions of the Pauli term, evaluated along the n~-lightlike direction, do 
not vanish completely when one employs the light-co ne gaug e — as opposed to the standard 



dx^A^ term. Moreover, it was shown in 



16 



24 



28 



-30 



35] that the transverse gauge field 



in the axial gauge at light-cone infinity is given by a total derivative, viz., 

1 f da + da~ e~ iq+ °°~~ iq ~ z+ 
A\oc~, z\ z ± ) = -^9j -j^y 2 ^ 2n5( q -) V <p(z ± ) , (20) 

whereas the longitudinal components are equal to zero. Thus, the field-strength tensor on 
the transverse segment vanishes: 

Fr(oo-,0+,£ ± ) = 0. (21) 

Therefore, expanding f[T6"j) . only the terms with longitudinal spin-dependent contributions 
survive, while those with J 7 (It) (or J 7 (la)) cancel out. Nevertheless, we verify the vanishing 
of these terms by explicit calculation in the next section. Hence, by virtue of Eq. ( 1201 ). 
expression ffTB"]) reduces to 



The bra-ket notation used will be explained in the next section. 
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poo poo 

[[oo*~, oo_l; oo~, Oj_]] • [[oo~, Ox; 0~ ; 0_l]].4+=o — 1 — W I dr I • Ailr) — ig drS'T(ur) 

Jo Jo 

POO pT 

-g 2 dr da (I- A(lr)) (I ■ A{la)) 
Jo Jo 

/*oo roo 

-g 2 dr do {I ■ A{lr))(S ■ 7{uo)) 
Jo Jo 

poo pr 

+ g 2 dr da{S ■ T{ut)){S ■ T{ua)) 



JO JO 

+ 0(g 3 ) . 

(22) 

IV. CALCULATION OF (VIRTUAL) GAUGE-FIELD CORRELATORS 

We are now able to calculate the spin-dependent contributions in Eq. ( fTUl) . which we will 
do up to the gf 2 -order level. Using light-cone coordinates (also in the transverse direction), 
the Pauli term reads 

S ■ T = S^ v (23) 

Imposing the light-cone gauge A + = 0, we obtain the following non-zero components of the 
field-strength tensor: 

= d+A- , F +l = d+A i , 
F~i = q-jv _ d i A - t jrH = d i A J _ d J A i _ (24) 

We proceed with the explicit calculation of the virtual gluon exchanges in Eq. ( jTOl) . rele- 
gating the inclusion of real-gluon contributions to Sec. I VI I The reason is that only the former 
are UV divergent and give rise to anomalous dimensions, while the latter contribute only 
UV-finite terms. To systematize the calculation of the various contributing correlators, we 
appeal to Table U in conjunction with Fig. [3j There are two different types of contributions: 
those terms in Eq. ( 1T61) which are proportional to g 2 stem from the evaluation of correlators 
between the standard gauge links and the enhanced ones. In Fig. [3] the latter are denoted by 
double lines with a ring attached to them in order to indicate the Pauli contribution which 
encodes spin effects. The standard gauge links are represented by simple double lines. The 

14 



other contributions to Eq. ( TT6l) . which are proportional to g, i.e., the terms Ux, M2, and U3 
in Table HI have to be contracted with the gauge fields generated by the Heisenberg fermion 
operators, cf. (j!7p . retaining again those terms which contribute to 0(g 2 ). It is understood 
that each of these terms has to be averaged over the fluctuations of the gauge field via a 
functional integration. This is done with the aid of Eq. (j27p using in what follows Dirac's 
bra-ket notation (...). 5 

The term (Ui) — graph (a) in Fig. [3] — reduces in the considered order of the coupling to 
what one obtains with the standard gauge links; it has been computed in our previous work 



in Ref. 24j. Term (U 2 ) — corresponding to graph (b) in the same figure — will be worked 
out below, whereas term {U3) vanishes by virtue of Eq. ( 12 ip . For the same reason, also the 
contributions termed (U 5 ), (U 6 ), (U 8 ), and (V(io) vanish as well. Moreover, it is proved in a 
few lines that (W 5 ) + (Uq) = 0. Term (U4) was computed in [24j and was found to vanish, 
while the term (W7), which represents the longitudinal selfenergy contribution of the Pauli 
term (graph (c) in Fig. [3]), will be computed further below; it amounts again to a vanishing 
contribution. Hence, the only remaining terms giving non-zero contributions are (U2) and 
(Ug). The first one stems from the interaction of the longitudinal gauge field, produced 
by the fermion, with the Pauli term along the enhanced longitudinal link — graph (b) in 
Fig. [3] — while the second one, represented by graph (d), describes the cross talk between 
the transverse gauge potential of the standard gauge link and the longitudinal part of the 
Pauli term (enhanced gauge link). Its calculation will be carried out below. Recall that 
the analogous cross talk between the longitudinal parts of the Pauli term and the standard 
gauge link vanishes because of Eq. (12TT) . 

Having sketched the general computational framework, let us now turn a spotlight on the 
calculation of the various terms, starting with (U4), while the Fermion- induced terms (Ui), 
(U2), and (Us) will be picked up in Sec. [V] The term {U4) represents the selfenergy of the 



24| . This can be seen from 



usual transverse gauge link and vanishes in the light-cone gauge 
the following equation 

^rt.^-X^^^l-AM^—, (25) 



Strictly speaking, one should write (...), 
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V 





V A 



(0 



(/) 



V 



V A 




FIG. 3: Feynman graphs contributing to the quark- in-a-quark TMD PDF. Double lines denote 
standard gauge links, while those supplemented with a ring represent enhanced gauge links with 
tensor-field (Pauli) contributions. Fermions and gluons are shown as solid and curly lines, respec- 
tively. Graphs (a), (b), (c), and (d) describe virtual gluon corrections; graphs (e), (f), and (g) 
represent real-gluon exchanges across the cut (vertical dashed line). The conjugated ("mirror") 
graphs are not shown for the sake of brevity. 



where we have used the Fourier transformation of the gauge field 

d"q 



(27T) 



(26) 



working in an u- dimensional momentum space (u = 4 — 2e). Employing this expression in 
Eq. (I25p and the gluon correlator in the light-cone gauge 

(27) 



(AMMq')) = H)C F (27r)V% + q>) D, v {q) 
in which the regularized free gluon propagator appears [cf. Eq. (jlj)], 



we get 



(Z4) = iC F l[ 



2nS(q-) 
q 2 + iO 



9^u 



dr \ do 
'o 



d 2 q ± e -^'^( r - CT ) 



2vr ./„ Jo J (2vr) 2 q\ + A 2 



(28) 



(29) 



One notes that the gluon propagator bears a pole-prescription dependence, codified by the 
symbol [q + ], whereas its non-zero parts are given by 

1 q i 



n = — D = i - 

g 2 + zO [q+] ' 13 q 2 + i0 



(30) 
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Despite the pole-prescription dependence of the gluon propagator, expression ( ]29l) does not 
depend on the pole prescription, because only the Feynman term of the gluon propagator 
contributes. Moreover, inspection of the last term in this equation reveals that it will be 
canceled by its mirror contribution anyway, i.e., finally, 



(W 4 > = 



(31) 



Note that this cancelation occurs in any case: polarized or unpolarized because there is no 
Dirac structure in this term. 

The next two terms (U5) and (Us), which contain expressions of the sort S ■ F{1±_t), can 
be treated in unison. To evaluate them we make use of the derivative of the transverse gauge 
field, viz., 



-iq + oo +iq±_ ij_(T lii 



(32) 



For the transverse gauge strength at light-cone infinity, one has 
implying for the Pauli term in the transverse direction 



(33) 



S ■ F{l ± a) = 2i 



(27T) 



-ig+oo +iqx-l_i_o i 

r e q± 



S +i A~(q) + S ij A 3 (q) 



(34) 



Then we find 



(W 5 > 



-2Cp— / dr I da 



T 



dq A 



d 2 q_i_ 



-iq_i_-l_j_(T-a) _ 



-s 



2tt J (2tt) 2 

+ S'iq' l{ 



An analogous calculation for the term (U 6 ) yields 

(U 6 ) = -(U 5 ) 



(35) 



(36) 



confirming that these two contributions cancel each other. 

Going forth, we can now compute the longitudinal self-energy spin-dependent (Pauli) 
contribution [graph (c) in Fig. [3] 



(Ur) 



dr / da {S ■ F(ut)){S ■ F(ua)) 
11 Jo 



(37) 
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using 



S . T {yTr) = 2i 



e~ iqUT q + 



S + - A-{q) + S- 1 A\q) 



(38) 



and employing the regularization embodied in Eq. (jlj) to obtain 

1 /.OO PT f ]UJ 1 

w - -4ic F - y * y d a y ^ 2 .%-) _ 

__2q- 



+ i0 



s + -s^ 



g—i g-\ — 2 5 j 



(39) 



It is easy to see that all three terms in the square bracket give vanishing results: The first 
term is zero because of S(q"). The second term gives also zero due to the oddness of the 
transverse integral, while the last one vanishes by virtue of 

1 



S~*S- = i ( 7 "7 1 ) 2 = 
[recall that = (1/4) [7^, 7^]]. Therefore, we finally get 



(U 7 ) = . 



(40) 



We consider now the term (14) in more detail and prove that it vanishes. As we already 
mentioned in connection with Table HJ this term, which represents the self-interaction of the 
transverse gauge links with the Pauli terms at light-cone infinity, vanishes by virtue of the 
particular form of the transverse gauge field in the light-cone gauge [see Eq. ( 120]) and the 
discussion below]. Here, we give a more detailed derivation of this result. By definition, this 
term reads 

/»00 PT 

(U 8 ) = dr da (S ■ T{It)) (S ■ JF{la)) . 
Jo Jo 

Therefore, we have 



(41) 



-4iC F - / dr da 

1 Jo Jo J (27T 



q 2 - A 2 + iO 



,■2? 



2q k 

[</+] " " ¥] 



S ik S n g ki 



(42) 



The term proportional to q~ vanishes by virtue of the delta-function. The second one is 
equal to zero because S q\q\ = 0. Taking into account that the Dirac structure of the last 
term can be rewritten as 



S tk S jk qi_q{ 



1 



1 



(43) 
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one obtains 



AtC F - I dr \ da [— f^e^H ql 



T 



dq + 







r — a 







(44) 



in agreement with the result presented in Table [H 

Consider next the mixed term (Ug), which expresses the correlation between the longitu- 
dinal Pauli term and the transverse gauge link [graph (d) in Fig. [3], viz., 



poo poo 

(Ug) = dr da (I ■ A(lr))(S ■ 7(ua)) . 
Jo Jo 

With the help of Eqs. (125)1 (|58"lh and (127)), we obtain 

1 fOO /-OO p Vj] 

' • 1 - ' y ^iq+oo - — iqj_-lj_cr— iq+ur 



(45) 



r*oo poo 

{Ug) = 2C Ff i 2 ^ I dr da 
1 Jo Jo 



x [S + -l l D^(q) + S-H i D i - 
which can be recast in the form 

"OO 



(27T) 



2tt 8{q~) g H 



(46) 



(Ug) = 2C F /i 



2c 



dr . . , 

T lo Jo J (2?r 



/ " e «'? + oo--i(j_ L -i_ L o--ic;+Mr 



2tt%2) 
g 2 + i0 



[S + -(q-l) + S-Tq+] 



(47) 



using Eq. fl30|) . Observe the important fact that the dependence on the pole prescription 
disappeared in the above equation on account of q + /[q + ] = 1, cf. Eq. ( I4"6"j) . As a result, this 
equation is valid for the advanced, retarded, and principal value prescriptions, as well as for 
the Mandelstam-Leibbrandt pole prescription, though it is not obvious that it holds true in 
general (see, e.g., Refs. 
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The r and a integrations in Eq. f T4T|) can be performed explicitly: 

dre~ iq+UT = — — — , 
o q + u — iO 

— i 

dae~ iq ^ 1 ^ = . 

q ■ I — iO 



(48) 
(49) 



Making use of the 
field at infinity 
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ojjowing relation, which stems from the structure of the transverse gauge 

(50) 



28 



29] 



- = 27ri S(q^ 

q + u — iO u 
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we get 

(Ug) = 2iC F fi 2€ — S + - / " . - n -= , (51) 

where the "gluon mass" A 2 was introduced in order to take care of infrared singularities in 
the gluon propagator. Taking into account that Tu = 1 and performing the q ± integral 

f d u ~ 2 q 1 _ % /4tt 

J (2tt)"- 2 q\ + A 2 - zO ~ 4^ V ^ 

we arrive at the following final result 

(W 9 > = -^C F [ 7 + ,7ir(e) ^Tr^y (53) 

that gives rise to a UV divergence. 

Its conjugated contribution, corresponding to the product of the gauge links 
[[oo~, £~, [[oo~, ooj_; oo~, amounts to the same expression (153"]) . i.e., 

(W 9 )t = (W 9 ) . (54) 

But there is a crucial difference: Now the Dirac matrix V in the definition of the TMD PDF 
stands on the right side of this expression — cf. Eq. ( 1T91) . Because the Dirac structure of 
Eq. ( |53l) is nontrivial, this will lead to different results. For instance, we get (using obvious 
abbreviations) 

(a) r unpo i. = 7 + : r unpo i.[7 + ,7~] = -[7 + ,7~]r u n P oi. , 

(6) r hclic . =7+ 7 5 : r helic .[ 7 + , 7 -] = -[ 7 + ,7ir hc i ic . , (55) 

(c) r trans . = icr 4+ 7 5 : r trans .[7 + ,7~] = -[7 + ,7 _ ]r t rans. 

From the set of these equations we conclude that after taking into account the conjugated 
(mirror) contributions, all the leading-twist two functions mutually cancel by virtue of the 
relation 

[ 7 + 7-] r tw _ 2 = -r tw _ 2 [ 7 +, 7 -] = 2r tw _ 2 . (56) 

This important property permits the probabilistic interpretation of twist-two TMD PDFs, 
because in every term 'tpTip, which behaves like a vector under ^-boosts, the pole contribution 
entailed by the correlation between the transverse gauge link and the Pauli term along the 
longitudinal direction disappears. 
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Remarkably, higher-twist distribution functions (e.g., twist three), behave differently, the 
reason being that they are characterized by a different Dirac structure that remains invariant 
under ^-boosts. For example, one has for T tw _ 3 = j l 

b + ,7l r tw _3 = r tw _ 3 [ 7 + ,7l , (57) 

so that the mutually conjugated contributions add to each other to give the net result 

r tw _ 3 (W 9 ) + (W 9 ) t r tw _ 3 = -g[ 7 + ,7ir(e) Uir^j (58) 

making it apparent that the pole contribution in that case is not vanishing. 

The last point that has to be verified is that the term (Wio) in Tabled vanishes. We shall 
do that without assuming the special form of the gauge field at infinity given by Eq. ( 1201 ) . 
Making use of the explicit form of the gluon propagator [cf. (|28|) ]. one obtains 

1 fV f, f d«q „ iq+(oo - u ^ iqi , lT 2ir6( q -) q + 



1 r°o poo /• juj 

(U 10 ) = -AiC F - J drj^ da J _|_ e ^(o°— )-^-^ 



x 



q 2 - X 2 + iO 

[q+] [q+] [q+] X 



(59) 



The first term equals zero due to the delta-function S(q~), while the third one vanishes by 
virtue of the antisymmetric-symmetric convolution S^q^q^ = 0. Performing the longitudi- 
nal line integral and taking into account Eq. fl50l) . that renders the last term vanishing as 
well, we reduce the above expression to 

w> = -uc^s->[ iT J ^ e «^_M_ . (60) 



Given that 



S+^qUi = ^ 7 + 7 - q\ , (61) 



we find after some standard calculations 



i poo p n 2 p iq±-lj_T 

W = V 7 - / dr 1 ' ^ 



T J J (2tt)-- 2 q\ + A 2 - iO 
= . (62) 

We thus conclude that the vanishing of the last term in Table [J can be proved even without 
additional constraints on the gauge field at light-cone infinity, like Eq. ( 1201 ). 
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V. FERMION (VIRTUAL) CONTRIBUTIONS 



We will consider now the terms (Ui) and (U2), pertaining to graphs (a) and (b) in Fig. [3j 
and also prove that the term (U 3 ) gives zero contribution. These are the O(g) terms in the 
expansion ( 12 2 p and have to be coupled to the fermion lines retaining their contributions up 
to the order g 2 . Term / has been considered in Ref. 24| and we will borrow the result from 
there: 

d"q 1 



(Wi) = C F 2m C^J 



(63) 



where the numerical factor C„ 



(2tt) w (p - q) 2 + iO q 2 - X 2 + iO ' 
{0;— 1;— 1/2} corresponds to different choices of the 



imposed pole-prescription in the light-cone gluon propagator (see Refs. [24 
singularity produced by this contribution, notably, 

1 



). The UV- 



uv 



— a s C F — i C Q 

e 



just cancels the prescription-dependent term in the UV-divergent part of the fermion self- 
energy graph in the light-cone gauge, bearing no relation to the spin-dependent part in 
question. 

The first novel contribution, ensuing from the Pauli term, is represented by the term (ZY 2 ) 
(graph (b) in Fig. [3]) making use of the notation we already employed in Eq. ( JT6l) . Hence 
the Pauli term with the tensor gauge field in the longitudinal direction becomes 

<Fq 



S ■ F{ut) 



2i 



(2k) 



-e - l i- UT q+ 



S + -A-(q) + S- l A l (q) 



(64) 



It stems from the interaction of the quark with the spin-dependent part of the longitudinal 
gauge link in the Pauli term. Consider first its longitudinal component, which we termed 
(Wf ), whereas for the transverse one we will use the notation Then, we have 



T 



dr 



d w q 
(2tt)^ 



e~ iq - UT 2K5{q-)q ] 



1 



1 



1 



(p-q) 2 + i0 q 2 -X 2 + iO 



(65) 



x [S + -(p - qh + 2q- + S-\p - g) 7 + gl] 
Taking into account that (7~7~ = , p* ± = 0) 



S~\p~q)l + q L ± = S- t [ 1 - 1 + (p + -q + )q i L + y 1 + qiqi] 



(66) 
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the pole-prescription-dependent term, containing cancels out and we get 

= C F7 - 7 + ^ fW [ f-^-^te-^- \ . (67) 

V 2 1 T J J (27r)"- 2 2vr p 2 - q\ + iO V ; 

Using the fact that uT = 1, we finally obtain 

K) = - 2 C^J W ±- 2W ^- . (68) 

Turning our attention to the conjugated contribution, we find out that the ordering of 
the Dirac matrices has changed: 

1 f°° r rl^n 1 1 1 

^ = ~ 2Cf t dT / 7^ ^ l9 UT ^(q-) q + - ^— - - — T 

-T Jo J (2vr) w (p - q) 2 + i0 q 2 - A 2 + 2O [g+J ( 69 ) 

x [-l + (p-q)S + -2q-- 1 +(p-q)S- l q i L \ . 
Therefore, one has 

,-U_ ^ -+--/■ ^ 2 «i 1 



Thus, combining these terms, we have for the leading-twist distribution the following 
final result 

(2^ 5( ^~ p2) 

= ^c F r tw _ 2 . (71) 

This contribution is UV finite and can be given a physical interpretation. Indeed, recalling 
Eq. ( fl6l) . we can express the above result in the form of a constant phase 

e~ iS « 1 - ^ 2 C F , (72) 

inherited to the TMD PDF by the Pauli term along the longitudinal gauge link. It is worth 
noting that this finding is valid not only for the unpolarized case with r tw _2 = 7 + , but 
also for the polarized case (r tw _2 = 7 + 7 5 or r tw -2 = iu t+ 'y 5 )- This is because the (leading) 
twist-two distribution functions are vectors under boosts along the ^-direction and their 
generic Dirac structure has the property 

r tw -27^7 + = 7 + 7~r tw _ 2 = 2r tw -2 • (73) 
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It is precisely this property that gave rise to the constant phase 5 and its cause can be traced 
to the correlation of the longitudinal gauge field, concomitant to the fermions, with the spin- 
dependent part of the longitudinal gauge link. This phase could, in principle, be absorbed 
into the soft factor R (cf. OH]))- However, this cannot be done in a universal way because 
the phase sign depends on the direction of the longitudinal gauge link, i.e., on the specific 
deformation of the integration contour via the it prescription. Inverting this direction, the 
phase factor f!72|) changes its sign and becomes 

e -iS _^ e+l s _ 

Therefore, the phases appearing in the SIDIS and the DY process turn out to have opposite 
signs: 

^sidis = -<W • (74) 

These UV features do not persist for the twist-three distributions. Indeed, their Dirac 
structures are invariant under z-boosts and behave like scalars, i.e., 

7"7 + r tw „ 3 =r tw _ 3 7~7 + , r tw _ 3 7 + 7~=7 + 7~r tw „ 3 . (75) 

As a result, the analogous expression to ( JTTT) now reads 

r t w-3<W 2 -> + (W 2 -)tr tw _ 3 = -^f j^[7 + ,7l ( 47r ^) r ( £ ) - ^} r tw-3 • (76) 

This quantity is UV divergent, meaning that the Pauli spin-dependent term will contribute 
to the anomalous dimension of the twist-three TMD PDF. 

We focus now on the interaction of the longitudinal spin- dependent gauge link and the 
transverse part of the gauge field originating from the fermions, namely, the term (ZY 2 ): 
I r°° r d^n 111 



V 27 T J J (2vr)- W ,H (p-q) 2 + i0 q 2 - \ 2 + iO 

Making use of the following simplifications of the terms with Dirac matrices, i.e., 



(77) 



s + (p - qh j q{ = s + ~ [7"(p + - q + ) + i + (p -q~) + i k Q k ] 

-+ -\h + ,rkl, (78) 
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and 



S-\p-q)Y = S~ l [r{p + -q + ) + l + {p--q-)+l k q k ]l l 

-> , (79) 

and proceeding along similar lines of thought as in the previous case, we find 

,±\ 1^ r..+ .-i fd^ 2 q± 1 



and 



so that 



W = -4*^17^73^ (80) 



r,,- 2 (z#) + (« 2 ± > , r„_ 2 = -jC F r lw ^ 2 . (82) 



The remarks which we have made in connection with (Z/f) apply equally well to Eq. (182j) . 
The computed phase is acquired through the interaction of the Pauli term along the longi- 
tudinal link with the transverse part of the gauge field accompanying the fermions and has 
to be added to the phase originating from the analogous interaction between the Pauli term 
and the longitudinal gauge field associated to the fermions. 6 

Thus, the full phase, ensuing from the interaction of the fermion fields with the spin- 
dependent (Pauli) term in the gauge links — crosstalk diagram (b) in Fig. |3] — is given accord- 
ing to (U2) by 

a 2 

e~ lS « 1 - % if- C F vr , (83) 

4:71 



where we have again taken into account Eq. (116]) . As we have already noted, this phase 
flips sign when the direction of the longitudinal link is reversed. As regards the twist-three 
distribution, we obtain 

rtw-a^) + (W 2 ± ) t r tw -3 = ^F^[7 + ,7l U-rr^) r(e)r tw „ 3 • (84) 

Comparison with Eq. (!76|) reveals that, taking their UV divergent parts together, their total 
contribution disappears leaving behind only a constant phase which, moreover, coincides 
with the one found for the twist-two distributions: <5 tw -2 = = ot s C-pii. 



6 The appearance of an imaginary contribution in the cusp anomalous dimension (unrelated to the Pauli 
term) was already discussed by Korchemsky and Radyushkin in [321 ] . 
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We complete our discussion of the fermion virtual contributions by considering the term 
(U3) in Table (TJ which describes the cross talk of the gauge field surrounding the fermions 
with the transverse spin-dependent gauge link. The discussion proceeds along similar lines 
as that of the previous term. Likewise, (U3) consists of two terms: (U^) and (W3 }. Consider 
first the contribution pertaining to the longitudinal gluons emanating from the quark fields: 



T 



dr 



d w q 
(2-nY 



-iq-oo +iqj_-lj_T 



(p-q) 2 + i0 q 2 -X 2 + iO 



[<?■ 



(85) 



This contribution vanishes because the first term equals zero by virtue of the delta-function 
5(q~), while the second one also reduces to zero due to the convolution S^q]_q]_. 

Continuing with the contribution from the transverse gluons produced by the quark field, 
we write 



-2CV 



1 

f 



dr 



d"q 



-iq-oo +iqj_-lj_T 



27i5(q- 



(p - q) 2 + iO q 2 - X 2 + iO 



S +i {p-q)l 



After performing the following transformations of the Dirac matrices 



S +l (p - g) 7 fc 5 +l 7 - 7 fc (p + - q + ) - S +l q k ± 



(86) 



(87) 



and making use of Eq. (12"0|) , one can recast Eq. 

l-e 



T \ iA-n 



00 

7 1 / dr 





into the form 

d u - 2 q 



(2vr 



7 



q ± — p 2 + iO 



e-l 




(88) 



Taking into account that ~ p + and T ~ p + , one sees that both terms in the square 
bracket are power suppressed. The first one, which is of 0{\p\/p + ) does not diverge; the 
second is of 0(p 2 /(p + ) 2 ) and is logarithmically divergent. In any case, both terms can be 
left out because we are only interested in the leading-twist contributions. Therefore, for our 
analysis the correlation between the transverse part of the Pauli term in the enhanced gauge 
link and the transverse gauge field produced by the fermion can be ignored. 
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TABLE II: Individual real-gluon contributions corresponding to the diagrams (e), (/), (g) in Fig. 



[3] and retaining terms up to order 0(g 2 ). 



Symbols 


Expressions 


Figure E] 


Uu 


J^drj^da (S-T(ut)) T (S - T(uv + £-,Z±)) 


(e) 


Ul2 


f™dr f °° da (I ■ A(lr)) r (5 • T(ua + f 5 


(f) 


Ul3 




(g) 



VI. REAL-GLUON CONTRIBUTIONS 



So far, we have presented the results of the calculation of the virtual gluon graphs, which 
(potentially) contribute to the UV-singularities of the TMD PDFs. Now let us turn to the 
real-gluon graphs, (e), (/), (g) in Fig. [3j The formal computation of these contributions is 
quite similar to that we already performed for the evaluation of the virtual graphs (b), (c), (d) 
in the same figure. 

The main differences are: 

(i) The discontinuity goes now across the gluon propagator, so that one has to replace it 
with the cut one. Then, in the light-cone gauge, we have 7 

~ q 2-\2 + t0 \ 9 + [g+] PV J 
DtocI^) = 2tt%+) 5(q 2 - A 2 ) (-gT + ^^^) • (89) 

(ii) The Dirac structures, abbreviated by T, stand now between Dirac matrices from the 
Pauli terms on different sides of the cut. 

(iii) The momentum delta-functions involve, apart from the "external" momenta p + and 
(k + = xp + , k±), also the "internal" loop momentum g M . 

The real-gluon contributions are listed in Table HT1 using analogous notations to those in 
Table U 

We start with the graph describing the interaction of two spin-dependent gauge links, as 



7 We omit here the discussion of the Mandelstam-Leibbrandt pole prescription 30, [37| , making the tacit 
assumption that the regularization of the [q + ] pole is q~ -independent. This will allow us to avoid an 
additional term in the cut propagator. 
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depicted in Fig. [3] (e): 



{Un) 



(2vr) 3 



e i(p+-fc+)r-i*x-€x / dr / dtT ( S . F(uT))r{S-T{ua+C;€x)) • (90) 



Making use of the cut propagator (1891 . and taking into account that the line integrals go 
along different paths so that they have not to be ordered, one has 



<Wii> 



T 



4:Cp— dr da 



d 4 q 



iiq+.U-iOy+iiq+u+iO)*^ _ x ) p +)£(2)( fc _ 



x(q + ) 2 2n5(q-) 6{q + )S(q 2 - A 2 ) 
2q- 



X 



s + - r ^ 



r s + - + 5+-r g ^ 



S-TS- j g ij 



(91) 



noting that the dimensional regularization becomes redundant in this case because all mo 
mentum integrals are finite due to the delta-functions. 
The d 4 q integration becomes trivial and hence we get 

k 



(s- 1 r s + - + s + - r s- 1 ) - — ±— - s~ts 

11 — x\p + 



We can now make use of the following formulas 



s~ l r s + - + s + - r s~* = -f 



S~TS~ l = -j 
2 ' 



(92) 



(93) 



(assuming, for instance, T = 7 + ) and perform an averaging with the help of ({TTj) to find 

Tr [{p + m){l + l5 s) 1 1 } = , 



Tr [(p + m) (1 + 75s) 7 ] = 2p 



(94) 



From this expression one may conclude that — within the given kinematics — the term (Z/q) 
is power suppressed. 

The next term represents the interaction between the longitudinal spin-dependent gauge 
link and the transverse gauge link at infinity — Fig. [3j graph (/). It reads 



00 ^00 



dr da{l-A{lT))T{S-T{ua + C^±)) 



(Wia) 



The further evaluation is analogous to that in Eq. (l4~6]) . giving the result 



(95) 



(Wi 2 > = C7 F 



T 



dr I da 



d 4 q 
(2^ 



5({1 - x)p + - q + )5 {2) (k x - q ± )2n 5(q~) q 



x e 



-iq+oo —iqj_-lj_cr—iq + UT 



2n5(q 2 - \ 2 )6(q + ) T [S+~l ■ q + S~ j r q + ] . 
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FIG. 4: Diagrams with a Pauli spin-dependent term, which may contribute to the time-reversal-odd 
phenomena. The intial/final-state boundary is denoted by a long vertical dashed line. 



The last term in the square bracket disappears due to the delta-function S(q + ), cf. Eq. (!50|) . 
Thus, one gets for the sum of diagram Fig. [3] (/) and its mirror counterpart a vanishing 
contribution: 

(u 12 ) + (w 12 ) t ~ r[ 7 + , 7 -] + [ 7 + ,7"]r = o . (97) 

We complete our task by calculating the contribution associated with diagram (g) in Fig. 
[3J which stems from the interaction of the Pauli term with the fermion field. Because it is 
of O(g), it has to be contracted with the gluon field in the quark-gluon interaction term 
if) Aip in order to contribute at order g 2 . The combined contributions graph (g) and its 
conjugate) of the longitudinal part are determined by the following combinations of Dirac 
matrices: 

(w 13 ) + (wr 3 ) t ~^r(p-g)7 + -7 + (p-g)r5- i . (98) 

Analogously, we find for the transverse part 

<*ti) + (W i L 3 ) t ~ S-r(p - g)Y - t* 03 - q)TS-* - S + -T(p - q)^ + ^ (p - q)TS + - . (99) 

After trivial manipulations with the Dirac matrices in the equations above, we finally arrive 
at 

(Kb) + (^r 3 ) t ~ ^ + (7~r + r 7 -) 7 l (ioo) 
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FIG. 5: The Feynman rules for the calculation of the one-gluon-exchange graphs shown in Fig. H]in 
the light-cone gauge and using enhanced gauge links that include the Pauli term. Rules are given 
for both sides of the final-state cut (long vertical dashed line). Vertical double lines represent the 
transverse gauge links at light-cone infinity, while the horizontal ones with arrowed rings are the 
spin-dependent light-like gauge links with the Pauli terms. 



and 



(z4) + (z4)t ~ -x P + [( 7 -r + r 7 -)y + 27T7-] 



101) 



From these results we conclude that these terms mutually cancel up to a power-suppressed 
correction. 

The main message from the computation of the real-gluon graphs containing spin- 
dependent terms is that they do not contribute to the TMD PDF in the leading-twist 
order. All physically important effects have their roots in the contributions of virtual-gluon 
exchanges. Exactly those diagrams are responsible for time-reversal-odd effects in more 
sophisticated models as we shall argue in the discussion to follow. 
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VII. SUMMARY AND CONCLUSIONS 



In this work we have presented a new gauge-invariant scheme for TMD PDFs which takes 
into account in the gauge links (Wilson lines) the Pauli term. This term describes the com- 
munication of spin degrees of freedom with the gauge field via the gauge-field strength and 
provides a rendering geared to calculational purposes of TMD PDFs for spinning partons. 
The key features of our approach can be summarized as follows: 

• The spin-dependent Pauli term, incorporated in the TMD PDFs as integral part of 
the gauge links, does not affect their UV-singular behavior in leading-twist order. 
Therefore, the structure of the gauge links in the soft factor, introduced before in 



24 



30| with the aim to define TMD PDFs in terms of gauge-invariant matrix elements 



with standard renormalization properties does not need to be changed. This proves 



the usefulness of the subtraction method, proposed in 
deal with rapidity divergences that cannot 
In fact, in our present analysis (and also in 
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34j , which provides a tool to 
)e controlled by dimensional regularization. 
we employed a soft renormalization 
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factor in order to cure overlapping UV and rapidity divergences and compensate this 
way the associated one-loop cusp anomalous dimension. 

However, the Pauli term contributes to the UV-divergences of the imaginary parts of 
the cut diagrams with virtual gluon exchanges. Though these effects cancel in the final 



result of the considered distribution of a quark in a quark, they signa 



ize that, within 



a more realistic context involving quark models with spectators |45|, 4f]], the spin- 
dependent terms may contribute to the interference diagrams, where these imaginary 
parts become crucial. 

By contrast, we found that the UV singularities of the higher-twist TMD PDFs (start- 
ing at twist three) are affected by the spin-dependent terms receiving contributions to 
their anomalous dimensions, which now become a matrix [see Eq. (|58|) ]. This is caused 
by an incomplete cancelation of UV divergences related to the fact that the z-boost 
induced by the Pauli term along the longitudinal link — pointing in one direction — is 
not counteracted by the conjugate contribution — pointing in the opposite direction. 
The net result is that only boosts and rotations around the transverse directions are 
left over and these give rise to a constant phase (see next item). Thus, to remedy the 
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definition of such TMD PDFs as densities, one has to compensate these divergences 
by introducing the Pauli term also into the soft factor. 

An important consequence of the presence of the Pauli term in the gauge links is 
that it gives rise to a phase entanglement, attributable to the interaction of this spin- 
dependent term with the companion gauge field of the fermion — diagram (b) in Fig. 
[3J 8 In technical jargon, the Pauli term along the longitudinal link generates z boosts, 
canceled by the conjugate link, and rotations along and around the transverse x and 
y directions, while the analogous term in the transverse gauge link produces boosts 
and rotations only along and around the transverse directions. The rotations of the 
longitudinal and the transverse link combine to produce a constant phase. It turns out 
that this phase correlation is the same for the leading twist-two and the subleading 
twist-three TMD PDFs, multiplying each of them as a whole. This means that ab- 
sorbing this phase into the soft factor for the leading distribution, the corresponding 
phase of the subleading functions is also removed, even though, as we explained in 
the previous item, these latter functions may lack a density interpretation. However, 
the Pauli-term- induced phase is not universal because it depends on the direction of 
the longitudinal gauge link. Reversing the direction of the gauge link, the phase flips 
its sign. Hence, it contributes with the opposite sign to the DY process relative to a 
SIDIS situation. This breakdown of universality indicates that the soft renormaliza- 
tion factor R does not fully decouple from the spin effects. For this to be the case, 
one would have to include into the definition of R spin-dependent terms (cf. Eq. Q)) 
and evaluate it along a topologically non-trivial contour (work in progress). 

To facilitate calculations with enhanced gauge links, we derive Feynman rules for both 
the left-hand side and the right-hand side of the final-state cut and display them 
in Fig. [5j The spin-dependent gauge-link propagator and vertices in the light-cone 
gauge are displayed in terms of double lines with arrowed rings around them. These 
feynman rules may be viewed as supplementing the set of Feynman rules given in 
4j for covariant gauges. Pay attention that the propagator of the transverse gauge 
link contains (in addition to the standard term originating from the line integration 



No such phase is induced by diagram (d), the reason being that this diagram does not involve a fermion 
line. 
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in momentum space) a numerical factor C^, which encodes the dependence on the 
pole-prescription — see Eq. (1201) . The cancelation of this dependence due to the soft 



factor is discussed in detail in our previous works 



30]. 



Let us now close our discussion by commenting upon possible consequences of the spin- 
dependent terms for models with spectators. The time-reversal-odd TMD PDFs, like the 
Sivers or the Boer-Mulders function, which are responsible for observable single-spin asym- 
metries (SSA)s, can be calculated by means of the graphs presented in Fig. H] — see for a 



recent analysis in 



47|. Such SSAs emerge as the result of the interference of the contri- 



butions of type (a) and (b) with their counterparts which bear no gluon exchanges. For 
instance, within the MIT bag model, non-vanishing time-reversal-odd TMD PDFs appear 
due to the interplay of the effects of the quark wave functions in the one-gluon-interference 
diagrams (see, e.g., Refs. 48-[si]). Therefore, in view of our results, one may conclude that 



the imaginary contributions (taken without their conjugated "mirror" counterparts), which 
derive from the spin-dependent gauge links, can affect the time-reversal-odd TMD PDFs 
even at the leading-twist level. For instance, the Sivers function of a quark having a flavor 
a is given by 



f^(x, k ± ) ~ \ff(x, k ± ; S ± ) - ff(x, k ± ; S A 



(102) 



so that it is defined by the sum of the imaginary parts of the diagrams (a) and (b) in Fig. 
HI This important finding and its phenomenological implications deserve further exploration 
and verification. 
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